BOUNDS ON COMPLEX EIGENVALUES AND 
RESONANCES 



A. A. ABRAMOV, A. ASLANYAN AND E.B. DAVIES 

Abstract. We obtain bounds on the complex eigenvalues of non- 
self-adjoint Schrodinger operators with complex potentials, and 
also on the complex resonances of self-adjoint Schrodinger oper- 
ators. Our bounds are compared with numerical results, and are 
seen to provide useful information. 



1. Introduction 

We consider self-adjoint and non-self-adjoint Schrodinger operators 
of the form 

Hf := -A/ + Vf 

acting in L 2 (R N ), where V is a potential vanishing sufficiently rapidly 
as \x\ — > oo. When V is a complex- valued potential we obtain bounds 
on the location of the eigenvalues of H, while in the case when V 
is real-valued we obtain bounds on its resonances. We also describe 
and implement methods for computing eigenvalues and resonances, and 
compare them both with each other and with others. 

There appears to be rather little literature on either of these prob- 
lems. The distribution of the eigenvalues of self-adjoint Schrodinger 
operators has been intensively studied for several decades, but most of 
the techniques are not applicable to non-self-adjoint operators. Sim- 
ilarly most of the literature concerning resonances is either computa- 
tional or considers the semi-classical limit and resonances close to the 
real axis. Our results are of quite a different character. 

The principal technique which we use to obtain bounds on eigenval- 
ues and resonances in sections 2 and 4 involves obtaining constraints 
on the numerical range of the operator H and of others closely related 
to it. In section 3, however, we describe a quite different technique 
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for ordinary differential operators which yields similar bounds but also 
enables us to prove that the number of eigenvalues of such operators 
is finite for potentials which decay rapidly enough. Section 5 is de- 
voted to obtaining a first order perturbation formula for resonances 
of a Schrodinger operator. For a general perturbation theory of reso- 
nances see, for instance, |I| and references therein. 

The remaining two sections concentrate on the computation of reso- 
nances. We have investigated such problems using a well-known com- 
plex scaling method for Schrodinger operators with dilation analytic 
potentials elsewhere 0, but here we consider Schrodinger operators 
in one dimension with rapidly decaying potentials which need not be 
analytic. In section 6 we describe two different methods of computing 
the resonances of such operators. In the final section 7 of the paper 
we apply our results to several examples. In particular, for comparison 
purposes we take operators for which the location of the resonances 
has already been determined. One of the examples of section 7 dates 
back to the series of papers |9|, [7|, ||, another to our earlier paper 0. 
We find that the two methods of section 6 have similar accuracy, but 
these methods are sometimes a lot much useful than the complex scal- 
ing method and sometimes considerably less accurate. In all cases the 
key numerical issue is the accurate computation of highly oscillatory 
integrals. 

2. Bounds on eigenvalues using complex scaling 

We initially assume for simplicity that V is a bounded complex- 
valued potential such that V(x) — > as \x\ — > oo. Since V is then 
a relatively compact perturbation of —A it follows routinely from || 
that the spectrum of H consists of R + := [0, oo) together with a finite 
or countable number of eigenvalues of finite multiplicity which can only 
accumulate at points of R + . 

We also assume that z — > V(zx) is complex analytic for all x G 
and all z in the closed sector 

{z G C : < 2arg(Y) < a} 

where < a < n and V{zx) — > as |x| — > uniformly with respect to 
z in compact subsets of the sector. We then put 

c{6) := ||V(e^ 2 .)||oc 

for all 9 satisfying < 9 < a. It follows from the three lines lemma 
that c(9) is a logarithmically convex function of 9. In the theorem 
below and elsewhere we adopt the notation 

B(a,r) := {z G C : \z — a\ < r}. 
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Theorem 1. Under the above assumptions on V the eigenvalues A of 
H which do not satisfy —a < arg(A) < all lie in the compact convex 
set 

fl {e- i0 R+ + B(0,c(9))} 

Q<6<a 

and their only possible accumulation point is at 0. 
Proof. If we put 

H z := -z~ 2 A + V(z-) 

then H z is unitarily equivalent to H for positive real z, so its eigenvalues 
are independent of z. The same applies by analytic continuation for 
complex z, in a sense which is familiar in the theory of resonances [[|. [| . 
When 

H e := -e- ld A + V(e i9/2 -) 

for < 9 < a then the essential spectrum of He depends on 9, being 
equal to e _ie R + . Since the eigenvalues can only accumulate in the 
essential spectrum and the only common point of the essential spectra 
of all of these operators is 0, this is the only possible limit point of the 
eigenvalues. 

In order to obtain bounds on the position of the eigenvalues, we need 
to introduce the numerical range 

Num(0) = {(#*/,/>: 11/11 = 1}, 

where we assume that all / lie in the independent domain of He. In 
the current situation Num(#) is a convex set which contains the entire 
spectrum of H e and is contained in 

e~ ie K + + 5(0, c(0)). 

It follows that if A is an eigenvalue of H then A lies in this set for all 
< 9 < a. □ 

Note If the sector of analyticity of V contains the positive real axis 
in its interior, then under similar assumptions one concludes that the 
only possible accumulation point of the set of all eigenvalues is 0. 

Other bounds of a similar type are available, depending on what 
information about the potential one wishes to use. 

Theorem 2. Under the above assumptions on V the eigenvalues A of 
H which satisfy Im(A) > all lie in the closed convex set 

P| {x + iy : x sin(#) + y cos(6 l ) < a(9)} , 

0<6Ka 

where the convex function a(-) on [0,a] is defined by 

a(9) := sup {lm(e i9 V(e id/2 v)) : v G R^} . 
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Proof. The main statement of the theorem depends upon the fact that 

Num(0) C {x + iy : x sin(0) + y cos(6>) < a(9)} . 

We now explain why a(-) is convex. If we define 

g(z,v) := lm{z 2 V(zv)} 

for < 2arg(^) < a, then as the real part of an analytic function 
g(-,v) is harmonic. Hence 

h(z) := sup{g(z,v) : v G R^} 

is subharmonic and constant on each ray through the origin. This 
forces a(-) to be convex. □ 

It is clear that similar bounds can be obtained for certain unbounded 
potentials. In the next theorem we state a general result of this type 
which uses quadratic form techniques. Here we make the same assump- 
tions on the analyticity of V in the sector as before, but replace the 
regularity assumptions as indicated. 

Theorem 3. Suppose that Vg := V(e %e / 2 ■) is relatively compact with 
respect to —A for all < 9 < a, and that for all such 9 the operator 
inequality 

\V e \ <le(-A)+(3 e 

holds where < 70 < 1 . Then every eigenvalue A of H which does not 
satisfy —a < arg(A) < lies in the closed convex set 

P| {C e + B(0,p e )} , 

0<6Ka 

where 

Cg := {z : — 9 — arcsin(76») < arg(^) < —9 + arcsin^)} . 



Proof. The only new component of the proof is the observation that if 
we put s := (—A/, /) then we obtain 

Num(^) C |J B(e- i6 s,sm( le )s + (3 e ) 
seR+ 

= C e + B(0,p e ). 

□ 
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3. Resolvent bounds on complex eigenvalues 

This section describes a completely different assumption under which 
one can obtain bounds on the complex eigenvalues, at least in one 
dimension. We replace the analyticity assumption on the potential V 
by the assumption that V lies in L 1 (R)flL 2 (R). This is known to imply 
that V is a relatively compact perturbation of —A acting in L 2 (R), so 
the spectrum of H = — A + V is again equal to R + together with 
eigenvalues which can only accumulate on the non-negative real axis. 
Moreover the domain of H is equal to W 2,2 (TV), which is contained in 
C (R). 

Theorem 4. Under the above assumptions every eigenvalue A of H 
which does not lie on the positive real axis satisfies 

\M < ll^ll?/ 4 - 



Proof. Let A = — z 2 be an eigenvalue of H where Re z > and let / be 
the corresponding eigenfunction, so that / G C (R). Then 

(-A + z 2 )f = -Vf (1) 

so 

-f = (-A + z 2 )- 1 Vf. 
Putting X := {V^ 2 , W := V/X and g := Wf G L 2 , we deduce that 

-g = W(-A + z 2 )- 1 Xg (2) 

so 

-1 G Spec(jy(-A + z 2 )-^). 

This approach is similar to the Birman-Schwinger principle for finding 
eigenvalues of self-adjoint operators (see, for instance, [[J). By eval- 
uating the Hilbert-Schmidt norm of the above operator, whose kernel 
is 

e -z\x-y\ 

W(x)^—X(y) 

we deduce that 

-2Kez\x-y\ 

\W{x) f \X{y)\ 2 dxdy>l. 

This implies ||1/|| 2 /4|z| 2 > 1. □ 

It appears more difficult to prove that the number of eigenvalues of 
H is finite or to obtain bounds on the number. The following theorem 
which treats the operator H in one dimension provides a result of this 
kind. 



6 A. A. ABRAMOV, A. ASLANYAN AND E.B. DAVIES 

Theorem 5. Let the potential V satisfy 

\\V{x)e rx \\ 1 < oo (3) 

for all 7 e R. Then the number of the eigenvalues of H is finite and 
all the eigenvalues (including real ones) satisfy 

\M<l\\v\\l 

Proof. Along with (j^) consider the equation 

u " - z 2 u = 0. 
Its linearly independent solutions are 

sinh(zx) 

Ui{X) = , U2{X) = cosh{zx). 

We shall seek the solution f(x) of ([I]) in the form 

f(x) = Ci(x)ui(x) + C 2 (x)u 2 (x), 
where the functions C\, C 2 satisfy 

C[(x) = V(x)f(x) cosh(zx), C' 2 (x) = -V(x)f(x f mh ( ZX \ 

We take C\{— oo) = z, C 2 (— oo) = 1 which corresponds to f(x) ~ 
'hen 

C x (x) = z+ ! V(t)f(t) cosh(zt) dt, 

J — oo 

c 2 (,) = i-f x vmm^lit. 



It follows from (^) that 

lim \Ck(x)\ < oo, k — 1,2. 



We define the function 

<p(z) =C 1 (oo) + zC 2 (oo). (4) 

As is easily seen, 

<p(z) = 2z + [ V(t)f(t)dt, (5) 

where 



R 



/» = /(x)e— = 1 + i f (1 - e- 2z M) V(t)/(t)d*. 

^ </-oo 

We observe that if A = — z 2 is an eigenvalue of H , then <£>(z) 
0, Kez > 0. For 2 in the right-hand half-plane we have 



h _ e ~ 2zr | < 2 Vr > 
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and, hence 



l/»l<i + n f l^)/»l d *- 

1*1 J — oo 

From this formula we get 





z\ 




3 






\V\\i 



sup \f(x)\ < 
This together with (|^) implies 

<p(z) = 2z + 1>(z), < 

Therefore <£>(z) 7^ for sufficiently large |z|. In particular, if \z\ > 

3 1 





z\ 


\\V\\i 


2 




-\\V\\i 



2 



V\\i then \ip(z)\ < 2|z| and cp(z) 7^ for all such z. 



Since /(x; z) is analytic in z, the function ip defined by (f|) is entire. 
Therefore the number of its zeros in any bounded domain is finite. The 
eigenvalues of H can only be found among the squares of the zeros of 
(p(z), hence the final result. □ 

Note Bounds similar to those in the above two theorems can be 
obtained for the corresponding operator on the half-line subject to 
the boundary condition a/(0) + bf'(0) = 0, \a\ + \b\ > 0, by minor 
adjustments to the proof. 

4. Complex scaling and resonance bounds 

We now apply the above results to the analysis of the complex res- 
onances of a self-adjoint Schrodinger operator H . We assume that V 
is real-valued, and that it can be analytically continued to the closed 
sector {z : —a < arg(z) < 0} with the same bounds on Vg as in theo- 
rem H itself has no complex eigenvalues but well-known arguments 
[|, |] show that those eigenvalues of Hg which lie within the sector 

{z : -0 < arg(^) < 0} 

are independent of 9. New eigenvalues may be emerge from the essential 
spectrum e~ l6, R + as 9 increases. The set of all such eigenvalues as 9 
varies over the permitted range is called the set of resonances of H. 
It is known that this definition is equivalent to others which do not 
depend upon complex scaling. 

We assume that Vg is bounded and vanishes at infinity for all 9 in 
the permitted range and put 

a(9) := sup {lm(e w V(e id/2 v)) : v G R^} 

as before. 

Theorem 6. Under the above assumptions those resonances of H which 
satisfy —a < arg(,z) < all lie within the closed convex set 

S :— O {x + iy : x sin(6 l ) + y cos(#) < a{9)}. 

0<9<a 
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Proof. Given 9, every resonance in the stated sector either has not yet 
been uncovered, in which case —a < arg(z) < —9, or it has been 
uncovered, in which case 

xsin(6 ) ) + ycos(9) < a(9). 

The proof is completed by a geometrical argument. The set 5* is the 
intersection of half-planes, and hence is closed and convex. □ 

The boundary of the set 5* is the envelope of the lines 

x sin(6 l ) + y cos(9) = a{9). 

which is given parametrically by 

x = a(9)sin(9) +a'(9)cos(9) 

y = a(9)cos(9)-a'(9)sin(9). 

We use this formula to find the form of the envelope as it emerges 
from the real axis. 

Theorem 7. // a(9) = a x 9 + a 2 9 2 + 0(9 3 ) for small 9>0 then a x > 
and a 2 > 0. Assuming a 2 > 0, the form of the envelope of S is 

V = -^r? 1 + o((x - a,f) 
4a 2 

for (x — ai) small and positive. 

Proof. An application of Taylor's theorem shows that 

ai = max{V(x) + x ■ VV{x)/2 : x G R^} 

and this implies that 

ai > M := max{V(x) : x G R^} > 0. 

The inequality a 2 > follows from the fact that the function a(-) is 
convex. The main statement of the theorem depends upon eliminating 
9 from the formulae 

x = a 1 + 2a 2 9 + 0(9 2 ) 

y = -a 2 9 2 + 0(9 3 ). (6) 

□ 

Note While one expects the hypothesis a 2 > of the above theorem 
to hold generically, the example of subsection |772] has a 2 = 0, and a 
higher order expansion would be needed to determine the form of the 
envelope near the real axis. 

Physical intuition suggests that there cannot be resonances A with 
Im(A) very small and Re A > M, which is confirmed by the example 
of subsection |0| (see figure 1), where Re A < M w 36.79. However, 
theorem |^ fails to prove quite as much. In the case in which V(x) = 
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x 2 e x l h for some b > 0, a simple computation shows that a\/M ~ 
1.256. 

A threshold for the resonances of an operator H is defined to be a 
real number 7 such that every resonance A of H satisfies Re A < 7. 

Corollary 8. Ifa>7c/2 then the resonances of the operator H possess 
a threshold. 

Proof. On choosing 9 = n/2, we see that every resonance satisfies 
ReA<a(7r/2). □ 

Corollary 9. Let r = 1 or 2 and let 

V{x) = [ e- s|x| "/i(ds) 
Jr+ 

for all x G R^, where /i is a finite signed measure on R + . Then 

1=1 H(da) 
Jr+ 

is a threshold for the resonances of the operator H. 

5. Perturbation formulae for resonances 

To obtain perturbation formulae let us assume that V(x) is real- 
valued and positive and rewrite the integral equation @ as 

-g(x)= W(x)G(x,y;z)W(y)g(y)dy, 



N 



JR 

where for all x, y the Green's function G(x, y; z) is analytic in z, Re z > 
0, positive for z > and satisfies G(x,y;z) = G(x,y;z). For certain 
values of N the kernel of (0) is known; the argument below does not 
use the explicit form of G. The above formula defines the operator 

A(z)g= I W(x)G(x,y;z)W(y)g(y)dy + g. 

JR N 

The resonances therefore are defined as such values of A = — z 2 that 
is an eigenvalue of A(z). The operator A(z) — I has a symmetric kernel 
for all z, Kez > 0, and is self-adjoint for positive real z. However, since 
we want to consider complex z let us decompose the operator A and 
the function g as 

A = A 1 + iA 2 , g = gi + ig 2i 

where the real and imaginary parts A\^ of A are self-adjoint, and g\$ 
are real-valued. The (non-linear in z) problem 

A(z)g = (7) 

is equivalent to the equation 

-A 2 At \f g 2 j _ Q (g) 



Ai A 2 I \ 01 
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The problem (§) involves a family of self-adjoint operators so we can 
now use standard perturbation theory techniques valid for such opera- 
tors. 

Suppose we have a resonance A = — z$, i.e. if we put z = z in (0) 
the problem has a solution g . Consider a perturbed operator A(z, e) = 
A(z) + eB(z) + 0(e 2 ) obtained, for instance, if we replace V by V e = 
V+eVi+0(e 2 ) in (JT|) (examples are to be found in section 7). Let us find 
the relevant resonance A = — z 2 assuming that z = zq+0(e). The latter 
assumption is justified by the above reasoning and by the fact that the 
operator function A(z) is analytic in z when z ^ 0. It is therefore 
expanded in powers of (z — Zq): A(z) = A(z ) + A'(zq)(z — zq) + . . 
the same is true for B(z). Following the conventional approach (cf. H]) 
we seek a solution of the problem Ag = in the form 

(A(z) + eB(z)){g + e gi ) = 0(e 2 ). 

Denote z — zq = ev + 0(e 2 ). Then the above equation can be rewritten 

as 

(A(z ) + eA'(z )u + eB{z ))(g + e gi ) = 0(e 2 ) 

yielding 

= CB(zo)go,3o) (g) 
(A'(z )g ,g ) { ) 

(here we use the property Ag = =3- A*g = 0). The operators 
involved are easily computed as 

A\z)g = f W(x) dG{x > V;z) W(y)g(y)dy, 
j r n az 

B(z)g= [ V 1 (x)G(x,y;z)W(y)g(y)dy + 

[ W(x)G(x,y;z)V 1 (y)9(y)dy. 

Finally, we formulate the main result of this section. 

Theorem 10. Let \q be a resonance of H = —A + V . Then the 
corresponding resonance of the perturbed operator H + eV% 

\ = \ -2v v ^>T e + O(e 2 ), 
v being given by ^) provided the denominator is non-zero. 

The first order correction to a given resonance is computed in terms 
of the solution of the corresponding unperturbed problem similarly to 
the well-known self-adjoint eigenvalue problem situation. Naturally, 
@ involves A' since the problem is non-linear with respect to z. 

The above approach essentially uses the symmetry of the kernel of A 
in the sense that (^) is reduced to the self-adjoint problem, which allows 
us to compute the correction v explicitly. However, a similar argument 
applies in a generic situation. Suppose, we have a non- self-adjoint 
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analytic family of operators A(z,e) whose simple eigenvalues X(z, e) 
are analytic functions of the two complex variables (see 0, Ch. VII, 
§4). To find the resonance of the operator A we put X(z, e) = 0; this 
equation can be resolved as z = z(e) provided || ^ by the complex 
analogue of the implicit function theorem. We therefore conclude that 
the resonance z depends analytically on the small parameter e. 

Perturbation formulae for a general analytic operator family are de- 
rived in |IJ, where the so-called resonance-eigenvalue connection is es- 
tablished in a very abstract context. In section 7 the author introduces 
an operator whose eigenvalues coincide with the resonances of the ini- 
tial operator; then in section 8 the perturbation formulae are obtained 
along the lines of |J. We, on the other hand, deal with the particular 
case of the Schrodinger operator, computing the first order term. The 
proof of theorem being fairly simple, permits one to switch to the 
self-adjoint problem, which is advantageous from the practical point of 
view when one wants to find the correction v. 



6. Rapidly decaying potentials 

We describe a technique for finding resonances which depends on 
rapid decay of the potential instead of analyticity assumptions. The 
method only works in one space dimension. We treat operators acting 
in L 2 (R + ) subject to generic boundary conditions at 0; similar methods 
apply on the whole line. 

6.1. Existence and uniqueness of the resonance solution. Reso- 
nances are conventionally defined to be solutions of the singular bound- 
ary value problem 

/" - (V(x) + z 2 )f = 0, x e R + (10) 
a/(0) + 6/'(0) = 0, |o| + |6|>0, (11) 
f[x) = e zx + o(\e- zx \), x->oo (12) 

Here z is a complex parameter and V(x) is a complex- valued function 
decaying rapidly enough at infinity: 

/ \V(x)\e 2zx dx < oo. (13) 

Lemma 11. For all z such that Kez > there exists exactly one so- 
lution of ( fIZ[ ) satisfying ftn\). 

Proof. Denote 

g{x) = f(x) - e zx , 

then 

g" - (V(x) + z 2 )g = V(x)e zx . 
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Let U-(x), u+(x) be two linearly independent solutions of (|l(J) defined 
by 

u±(x) ~ d±e ±zx , x — ► oo, (14) 

d± being some constants. 

We look for the function g(x) in the form 

g(x) = C-(x)u-(x) + C + (x)u+(x), 

where the functions C_, C + solve 

C'_u_ + C>+ = 0, 

C'_u'_ + C>; = V(x)e zx . 

In other words, 

1 r°° 

C ± (x) = T W J V(t)e zt u T (t)dt. (15) 
Here the Wronskian W = 



U- u. 



v! u'. 



and (O) guarantee the convergence o: 



is constant. Assumptions ( [131) 
the integrals in flT5| ) and imply 



\C + (x)\ = o(e- 2zx ), |C_(x)| = o(l), x -> oo. 
This, in turn, means that 

= o(e~ zx ), x^oo (16) 

as required. The way we construct this solution also implies its unique- 
ness. □ 

Now that the lemma is proved we are led to the following problem: 
find such values of z that the equation ( 110 ) has a solution satisfying 
both of the boundary conditions flTT| ) and (|P2|). The values of A = —z 2 
are the resonances of our problem. 

6.2. First method. Given any z we can compute at x = the value 
of g{x;z) satisfying (|l6|) . The resonances are determined by 

a(l + g(0;z)) + b(z + g'(0;z)) = 0. (17) 

The resonance solution f(x;z) then satisfies the boundary condition 
©■ 

One possible method of evaluating g(0), g'(0) is as follows. Denote 

a±(x) = u'±(x)/u±(x) =F z- 

These two functions can be easily computed as the solutions of the 
Cauchy problems 

a'± + a4 ± 2a±z - V = 0, (18) 

tt _(oo) = 0, a+(0)=a . (19) 
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Here the solutions of the problems (0), Q19D are assumed to be bounded 
on (0, oo). We then have 



o 



U-(x) = m_(0) exp ( / a_(t)dt — zx 

and, hence 

C+(x) = -W _1 «_(0) exp ( / a_(t)dt ) J+fc), (20) 

where 



I Ax) 



Similarly, 



where 



/OO / POO \ 

V(t)exp(-J «_(r)drjdt. 



C_(x) = W _1 u + (0) J -(^) ( 21 ) 



"OO / 

J_(x) = / V(t) exp ( 2zt + / a + (r)dr ) dt. 



.r 



A simple calculation gives us 

W = m_(0)m + (0)(« + (0) - a_(0) + 2z). 



We can now evaluate g(0) and </(0) using (P0|) , (|2l|). The final 
formulae are 

ro s = /-(0)-/+(0)exp(/ R+ a_(x)dx) 

= /-(0)(a-(0) - z) - I+(0)K + g) exp(/ R+ q_(a:)dx) 
^ J a_(0) + 2z ' 1 j 



o 



We did not specify the initial value aq_(0) in fll8|) ; now it is clear that 
we can take any a 7^ a_(0) — 2z. Computing g(Q), g'(0) by formulae 
(p2|), fl23|) for a range of 2 we then solve the equation (|T7|) by a standard 
iterative procedure and, finally find the resonances. 

Following this procedure we have to integrate flT8|), ( |IT?D numerically. 
The equation for a_ is solved from right to left, for a + from left to right, 
which means the integration is stable in both cases. As we compute 
a±, at each integration step the integrals J* a + (t)dt, J x °° a_(i)d£ and 
I-(x), I+( x ) are successively evaluated. Note that the integral J_ may 
converge rather slowly depending on the rate of decay of the potential 
V(x) and on the size of \z\. One has to choose a suitable method for 
evaluating J_ according to the behaviour of the integrand. 
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6.3. Second method. Let us normalise the formerly introduced u-{x) 
so that u^(x)e zx — ► 1 as x — > oo (gL = 1 in Then 



u_(x) = e~ zx ^(x), 7(x)=exp^y a-(t)dt J . (24) 
A general solution of ( [TOD is given by 

f(x)=u_(x)(^c l + c 2 ^ty j » ( 25 ) 

where Ci, C2 are arbitrary constants. We take C2 = 2z and look for such 
Ci that /(#) = e 2a; + o(w_(x)), x — > oo. From (|25|) we then have 

dt 



Ci = lim — — — 2z 



U-{x) Jq U-(t) 



\2 



= B _(o)- + / °"(«-W-« , -W)-^ dj = 

T e^( 7 (*)(2,- «(»))- 2^ 

Everything is expressed in terms of a_ in the above formula; u_(0), j(x) 
are evaluated by means of (^). Having found ci we observe that 
/(0) = ci«_(0), f (0) = cim'_(0) + 2zu-(0)-\ It remains to find the 
resonances from (|11|) as before. 

The two methods described above along with the method of |2j based 
on complex scaling have been successfully applied to several problems 
to be considered in the next section. For a number of problems dealt 
with here (see the next section) the numerical efficiency of the two 
methods appears to be approximately the same. Choosing between 
them is a matter of convenience, although there may be situations 
when one of them is somewhat better than the other. We observe that 
the methods of this section are equivalent within the range of their 
consistency. 

7. Examples 

7.1. Gaussian potential. In certain cases no useful information about 
the resonances of the operator H can be obtained by the complex scal- 
ing method. For example if we put 

H g f(x):=-Af(x)-e-* 2 f(x) (26) 

acting in L 2 (R + ) subject to Dirichlet or Neumann boundary conditions 
at 0, then complex scaling can only be applied for angles 9 < ir/2. 
However all known resonances A of this operator have —n < arg(A) < 
— 7r/2, so they cannot be computed using the complex scaling technique 



We take the operator H g as an example and find its resonances by 
means of the method of subsection 6.2. The results given in table [I] 
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Table 1. Resonances for the Gaussian potential 



n 




1 


(D) 


-0.52964412 


2 


(N) 


— 1.23692584 - 


- 3.48426766? 


3 


(D) 


— 1.33981054 - 


- 6.70063305? 


4 


(N) 


— 1.45338641 - 


- y.yt>21t>59Z2 


5 


(D) 


— i.jy^uyiou — 


1 'i 1 zLQRQzLQS? 
10. i^tcjuy^zofc 





f'lVA 

(IN) 


-1.69788643 - 


16.29120850i 


7 


(D) 


-1.75150063 - 


19.43849866i 


8 


(N) 


-1.80493057 - 


22.61135566i 


9 


(D) 


-1.87573378 - 


25.77287400i 


10 


(N) 


-1.93288884 - 


28.9110644H 


11 


(D) 


-1.96680887 - 


32.05477086i 


12 


(N) 


-2.00314399 - 


35.21335337i 


13 


(D) 


-2.05016204 - 


38.36530926i 


14 


(N) 


-2.13212307 - 


41.49620539i 


15 


(D) 


-2.11478618 - 


44.64676494i 


16 


(N) 


-2.14162017 - 


47.79467138i 


17 


(D) 


-2.16811153 - 


50.85272557i 



are obtained up to the accuracy of 10~ 8 . Some of the smaller reso- 
nances of the same operator, presented very approximately in [TIJ , are 
at least in a qualitative agreement with our results. The labels (D) and 
(N) in table p] relate to Dirichlet and Neumann boundary conditions 
respectively. 

Note The operator H g is also known to have real eigenvalues. A sim- 
ple variational computation shows that there is at least one eigenvalue 
A ~ —0.335 corresponding to the Neumann boundary condition at 
x = 0. On the other hand, it follows from JD], p. 91] that H g cannot 
have more than one negative eigenvalue. This implies that the first 
entry in table [I] is indeed a resonance with a negligibly small imaginary 
part. 

As well as the Gaussian potential we consider its perturbations 

V £ (x) = - exp (e-\l - (1 + 2EX 2 ) 1 ' 2 )) , 

for small e > 0. The reason for making this rather complicated choice 
of perturbation is that although one cannot use complex scaling for 
the original potential for angles bigger than tt/2, the approximating 
potentials allow one to use complex scaling with angles up to tt, which 
is enough to capture the resonances of this operator. The convergence 
of the first two resonances of the perturbed operator to those of H g as 
e — > is shown in table 0. The table provides a numerical confirmation 



of the result of theorem |10[ The resonances of the perturbed and the 
original operator agree with the perturbation formulae of section [5]. 
The difference is seen to be of order 0(e) as expected. 
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Table 2. Perturbed resonances 



£ 


- A« 


Af - A( 2 ) 


5 • 

1 • 1(T 4 
5 • 10- 4 
1 • 10" 3 


1 • 10" 4 
2 • 10~ 4 

i • icr 3 

2 ■ 10" 3 


-2 • 10" 4 + 6 ■ 10" 4 i 
-3 • 10~ 4 + 1.1 • 10 3 i 
-1.6-10- 3 + 5.3-10" 3 i 
-3.2 • 10~ 3 + 1.06 • 10" 2 i 



7.2. A slowly decaying modified Gaussian potential. The oper- 
ator 

H b f{x) := -Af(x) + x 2 e-* 2 / b2 f(x), (27) 

b being a real constant, has been studied in [Q. Here we provide more 
detailed information on its resonances for b = 10. Corollary 8 is not 
applicable to the operator Hi, but it seems from the numerical evi- 
dence that a threshold exists. The resonances computed numerically 
are shown in figure 1. We used the methods of section 6 to find some 
of them and compared the results with those obtained in || by means 
of complex scaling. The methods agree with each other for Rez < 40; 
for higher resonances complex scaling appears to be the best out of the 
three methods. 

Theorem ^| can be applied to obtain a bound on all resonances within 
the sector 

{z : — 7r/2 < arg(z) < 0}. 

For this example we have 

a{9) = b 2 sup{se- scos(9) sin(2fl - s sin 9) : s G R + } 

for every < 9 < ir/2. This may be used to compute the envelope of 
the set S using the parametric formulae (§); this is compared to the 
actual resonances (found by complex scaling) in figure 1. Of course a{9) 
may be evaluated purely numerically in more general cases. Naturally, 
the bounds on the resonances do not compare closely with their actual 
position, but it should be noted that the bounds are obtained with 
much less computational effort and do provide useful qualitative and 
quantitative information. 

7.3. The potential of Rittby et al. We finally study the potential 

V(x) = (x 2 - J)e- 0Ax2 + J, J = 1.6. 

In H [7|, |] one can find controversial remarks on the related resonance 
problem and two quite different series of resonances for the same op- 
erator. Our methods give broadly the same resonances as presented in 
It appears that the boundary condition at infinity taken in [0] is 
different from ( |T2"D as has been mentioned in |3J, hence the difference 
between the results. 
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Figure 1. Resonances of H b} b = 10, and the set S 

By means of the methods of section |6| we have computed about 20 
resonances, which coincide with those presented in up to 4 decimal 
places. Computing higher resonances by either of the methods of sec- 
tion ^| turns out to be quite difficult, whereas the method of complex 
scaling still works well. Indeed, using the technique of 0] we obtain 
stable results which agree very well with all the resonances quoted in 
[[]], whose method is similar. This suggests that if a potential decays 
fairly slowly and the complex scaling method is applicable, then it may 
often give more accurate results than the methods of section 6, es- 
pecially for higher resonances. We have already commented that for 
rapidly decaying potentials like Gaussian, the methods of section || can 
be efficiently applied. It is possible to improve those methods by using 
an advanced technique for computing the highly oscillatory integrals 
involved (see the concluding remark of subsection |6]^) . 

Apart from the resonances found in [^, we have also discovered 6 
more resonances: two near each of the points Ai = 0.69 — 7.91i, A2 = 
1.26 — 8.51i, A3 = 2.08 — 11.61i. Each pair includes an even and an odd 
resonance very close to one another. This closeness hinders accurate 
computations and is probably the reason why those resonances are 
missing in |], |J. None of our three methods allows us to locate them 
accurately, but there is clear evidence that they exist. 
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